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Abstract. The spin-down mechanism of accreting neutron stars is discussed with an application to 
one of the best studied X-ray pulsars GX 301-2. We show that the maximum possible spin-down 

torque applied to a neutron star from the accretion flow can be evaluated as — jX 2 / (r m r CO r) 3 ' /2 - 
The spin-down rate of the neutron star in GX 301-2 can be explained provided the magnetospheric 
| radius of the neutron star is smaller than its canonical value. We calculate the magnetospheric radius 

considering the mass-transfer in the binary system in the frame of the magnetic accretion scenario 
suggested by V.F. Shvartsman. The spin-down rate of the neutron star expected within this approach 
i-C ' is in a good agreement with that derived from observations of GX 301-2. 
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INTRODUCTION 

Observations of the long-period (P s ~ 685 s) X-ray pulsar GX 301-2 have shown periods 
(lasting up to a few years) when the neutron star spin frequency (v = 1/P S ) is steadily 
decreasing at the rate Vo ~ — 10~ 13 Hzs _1 [1]. This pulsar has been identified with a 
neutron star in a High Mass X-ray Binary with the orbital period P OI \, ~ 41.5 d. The 
\ massive component is the early-type supergiant (Wray977 D2|]). It underfills its Roche 

lobe and looses material at the rate M Q ut — lO^MQyr -1 in a form of relatively slow, 
v w ~ 300 — 400km s _1 , stellar wind dH]. The neutron star moving in the wind of its 
companion captures material and accretes onto its surface at the rate: 

>< : M a = ^^~5xl0 16 L 37j R 6 m- 1 gs- 1 , (1) 

where — R ns / 10 6 cm and m = M ns / 1 .4Mq are the radius and mass of the neutron star, 
andLr? =L X /10 37 ergs 1 is the X-ray luminosity of the pulsar normalized according 
to |4j,[3|]. Observations have revealed no signs of developed Keplerian accretion disk in 
the system. Besides, the relative velocity of the neutron star in the wind of its massive 
companion, v re i = v w + v ns , exceeds the upper limit (see D50), 

v cr * 200 ^ 2 4 M3o 1/14 m 11 / 28 MT (^) ^ kms- 1 , (2) 

at which the angular momentum of the captured material is sufficient to form a Kep- 
lerian disk. Therefore, the accretion process in GX 301-2 is usually treated in quasi- 



spherical approximation. Here v ns ~ 250km s _1 is the linear orbital velocity of the neu- 
tron star, /I30 is its dipole magnetic moment in units of 10 30 Gcm 3 , A/17 is the rate at 
which the neutron star captures material from the stellar wind in units of 10 17 gs _1 and 
<d0.2 = <d/0-2 is a parameter accounting for the angular momentum dissipation due to 
inhomogeneities in the accretion flow. It is normalized to its average value according 
to results of numerical simulations of wind-fed accretion in the approximation of non- 
magnetized accretion flow (see [6] and references therein). 

Interpretation of the observed spin-down of the neutron star in the frame of popu- 
lar scenarios of spherical accretion encounters some difficulties. Lipunov (see |7|, |8Q) 
has shown that the spin-down torque applied to the neutron star undergoing quasi- 
spherical accretion can be evaluated from the expression 1C$ = fc t jU 2 /r 3 or . Here r cor = 

1 /3 

(GM ns / Gt) 2 ) is the corotation radius of the neutron star spinning at angular velocity 

(O s = 2n/P s , and kt < 1 is the efficiency parameter. The spin-down rate of the neutron 

(si 

star in GX 301-2 in this case is limited as 27r/v s ^ sd and its absolute value, 

v s < 3 x 10- 16 k t I 45 l m- l R 6 6 P 6S 2 5 Hzs- 1 , (3) 

turns out to be a factor of 340 less than that derived from observations. Here ^45 = 
7/10 45 gcm 2 is the moment of inertia of the neutron star, i>685 is its spin period in the 
units 685 s and #crsf = 4 x 10 12 G is the magnetic field strength on the stellar surface 
estimated through observations of the cyclotron line in the X-ray spectrum of the pulsar 

a. 

The value of spin-down torque in the scenarios accounting for an influence of the X- 
ray luminosity of the pulsar on the accretion flow [10] and a possibility for a hot gaseous 
envelope to form around the neutron star magnetosphere [jll| , under certain conditions 

can be as large as = zoMcCO^. Here M c — nR^PooV^i is the mass of gas with which 
the neutron star interacts in a unit time moving through the wind of its companion, 
Rq = 2GM ns /V r 2 1 is the gravitational capture (Bondi) radius and is the density of 
stellar wind at the Bondi radius. r m is the radius of the neutron star magnetosphere 



which in the case of spherical accretion is r m > r a , where IU2I1 

2 \ 2 / 7 



M c (2GM ns ) 1 /2 



(4) 



and zq < 1 is a numerical parameter. The expected spin-down rate of the neutron star in 
GX 301-2 in the frame of these scenarios 

/ n \ 8 / 7 

v«= ! 2xlO-»Hz S -' z „/ 45 '^m- 5 /^f 7 Z.^(-^) , (5) 

\#CRSF/ 

exceeds an estimate derived from the canonical model of quasi-spherical accretion (see 
expression |3}, but still remains a factor of 5 less than the observed value. 

Analyzing this discrepancy, Doroshenko et al. [1] have assumed that the magnetic 
field strength on the surface of the neutron star in GX 301-2 is a factor of 30 higher than 



that estimated through observations of the cyclotron line in the X-ray spectrum of the 
source, #crsf- The neutron star under this assumption proves to be an accreting mag- 
netar with the surface field in excess of 10 14 G. In this paper we argue that so dramatic 
revision of our notions of magnetic fields in accreting neutron stars is premature. An 
apparent contradiction between the predicted and observed spin-down rate of the neu- 
tron star indicates that the accretion picture in the long-period X-ray pulsars may differ 
from commonly adopted scenarios. In particular, the observed spin-down rate of the 
neutron star can be explained in the case of magnetic accretion [| 1 3|1 without additional 
assumptions about the value of the neutron star magnetic field. 



MAGNETIC ACCRETION 



As was first shown by Shvartsman [|13ll . a picture of the wind-fed accretion in a massive 
binary strongly depends on the value of magnetic field in the matter captured by the 
compact star. The magnetic field in the free-falling accretion flow is dominated by 
the radial component. This is connected with the fact that in the process of spherical 
accretion the transverse scales contract as ~ r~ 2 , while the radial scales expand as 



r V 2 [14]. Therefore, under the condition of the magnetic flux conservation the field 
strength in the accretion flow can be approximated as B T r~ 2 fl 1 311 . This implies that 



the magnetic energy density, E m = B^/Sn, in the free-falling material increases as it 
approaches a compact star as 

E m {r)=E m {R G )(—) , (6) 



r 



while the kinetic energy of the accretion flow is 



(R \ 5/2 

£ram(r)=£ram(*G) — • (7) 



V r 

Here E m (Rc) = P^E^Rq) is the magnetic energy density in the accretion flow at the 
Bondi radius normalized by the thermal energy of the flow, E lh (R G ) = p^c^Rq), with 
the parameter /3, and c s (Rq) is the sound speed in the captured material. E mm (RQ) — 
pooV 2 el is the kinetic energy density in the flow at the Bondi radius. 

Then solving the equation E m (R s h) = E Tam (R s \ i ) we can find the distance 7? s h (here- 
after, Shvartsman radius) at which the magnetic energy density in the accretion flow 
reaches its kinetic energy density in the form 

Shvartsman [fl3ll has pointed out that for the accretion flow in the spatial region r < R s ^ 
to proceed further the dissipation of its magnetic field should occur. Otherwise, the 
magnetic energy in the flow would exceed its gravitational energy which contradicts 



the energy conservation law. Thus, the radial velocity of the accreting material in this 
region can be expressed as v r ~ r/t rec , where 



r\mVA \V A J 

is the characteristic time of field dissipation due to magnetic reconnection. Here Va = 
B r / (Aitpy/ 1 is the Alfven velocity in the accretion flow, tff = r/Vg and Vjj(r) = 

(GM ns /r)^ 2 is the free-fall time and velocity. r\ m is the efficiency parameter of the 
reconnection process. Its value depends on the physical parameters of plasma and 
magnetic field configuration in the region of reconnection and ranges in the interval 
< r\ m < 0.1 fl^l . Since Va < Vff (the equality is met at the Shvartsman radius), the 
timescale of field dissipation in the flow remains significantly less than the dynamical 
(free-fall) time, %, during the entire accretion process. This confirms the validity of as- 
sumption about the magnetic flux conservation in the free-falling material, but, on the 
other hand, this implies that the flow is decelerated by its own magnetic field at the 
Shvartsman radius and accretion switches into the diffusion regime. 

Basic conclusions of the scenario suggested by Shvartsman have been later confirmed 
with quantitative assessments by Bisnovatyi-Kogan and Ruzmaikin lUl 0], and by 
results of numerical simulations of the spherical magnetic accretion presented in ||19[ 
20Tl . These authors have shown that magnetic field amplification in the free-falling 
material leads to deceleration of the accretion flow at ~ R s ^ and its shock-heating up to 
adiabatic temperature. Accretion inside the Shvartsman radius occurs on the timescale 
of magnetic field dissipation. However, the accretion picture essentially depends on 
the efficiency of cooling processes in the accreting material. If the cooling time at 
the Shvartsman radius, t coo \(R s h), exceeds the heating time due to magnetic energy 
dissipation, t KC (R s h), accretion inside the Shvartsman radius occurs from a hot turbulent 
envelope with some portion of matter leaving the system in a form of jets [fl9i l20ll . 
Otherwise, the flow is transformed into the magnetic slab with parameters depending, 
in p articular, on the geometry of the large-scale magnetic field in the accretion flow 



Vff 



in p arti 



MAGNETIC ACCRETION IN X-RAY PULSARS 

A possibility for magnetic accretion to be realized in an X-ray pulsar is determined by 
relation between the Shvartsman radius and the radius of the neutron star magneto sphere. 
If R s h < r a , the accretion flow reaches the magneto spheric boundary of the neutron star in 
the free-fall regime and the accretion process in this case can be considered in the quasi- 
spherical approximation. Otherwise, the accretion picture should be treated in terms of 
magnetic accretion. Solving inequality i? s h > r a for v re i, we find that magnetic accretion 
in the massive X-ray binaries can occur provided v re i < v mca , where 

v mca = /3- 1 / 5 (2GM ns ) 12 / 3 V 6/35 M C 3 / 35 c 2 J 5 . (10) 



For typical parameters of the long-period pulsars (see [0]), this yields 

c 680 />-"* m ' 2/ V„ 6/35 Mf 7 /35 do^) 2 ' 5 (1 1) 

Thus, for moderate values of the parameter /3, v mca significantly exceeds v cr (see Eq.©. 
This allows us to distinguish a subclass of X-ray pulsars in which the magnetic accretion 
scenario can be realized. This subclass is determined by the following condition v cr < 

^rel ^ ^mca- 

The parameter /3 in the material captured by the neutron star at the Bondi radius can 
be, in the first approximation, estimated as follows. Let 5 ms denote the magnetic field 
strength on the surface of the massive counterpart to the neutron star. The field strength 
in the stellar wind ejected by this star is decreasing as it is moving away from the star 
as B w (a) ~ B ms (R ms /a) , up to the distance at which the kinetic energy density 
of the wind becomes equal to the magnetic energy density of the stellar dipole field. 
Here R m . is the radius of the massive star, and a is the distance from its center. The 



magnetic field strength in the stellar wind in the region a > a\ decreases cT 1 plQ . 
Thus, the magnetic energy density in the stellar wind in the vicinity of the neutron star 
can be evaluated using the expression E m (ao) = jU^/ (ina^afy, where ao is the orbital 
separation. Putting the numerical values yields 

^-O.SSergcm-^r^^^)^^-) , (12) 

where = ao/10 13 cm, and /i ms is the dipole magnetic moment of the massive com- 
panion which was normalized according to spectropolarimetric data on massive stars 
11221 123U24I1 . These observations have shown that magnetic field strength of a few thou- 
sand Gauss is not unique among early-type stars. 

The thermal energy density in the material captured by the neutron star at the Bondi 
radius, E^(Rg) = P<x>c 2 , can be expressed by taking into account that = M c / TlR^v K \ 
as follows, 

E th ~ 0.02ergcm~ 3 mT 2 M X1 ( — — T ) ( — =Y (13) 

th e 17 \500kms- l J VlOkms-V V ' 

These estimates demonstrate that the parameter /3 in the stellar wind captured by the 
neutron star can be of the order of or even less that unity in not too wide pairs with 
a moderate stellar wind velocity. It should be noted that /3 ~ 1 has been measured in 
the solar wind plasma in the vicinity of the Earth. This indirectly confirms the validity 
of assumption about significant magnetization of stellar wind adopted within magnetic 
accretion scenario. 

Finally, transformation of the spherical flow to the magnetic slab inside the Shvarts- 
man radius 11171 Il8ll can occur provided t coo \ < t KC . The most effective cooling mech- 
anism in the accretion flow in X-ray pulsars is the inverse Compton scattering X-ray 
photons emitted from the neutron star surface on the hot electrons of plasma surround- 



ing its magnetosphere [1121 125|]. Solving inequality t c (R s h) < f rec , where 



. , 37tr 2 m e c z 

•< (r) = -r^r <14) 



is the Compton cooling time [25], we come to the conclusion that a formation of the 
magnetic slab surrounding the neutron star magnetosphere is possible only if the X-ray 
luminosity of the pulsar exceeds the critical value: Lx > L cr , where 

icr = 3 x 1033 „./4 R -m ( ^ 1/2 ergs ->. (15) 

Here m e is the electron mass, and Oj is the Thomson cross-section. 

The plasma density at the inner radius of the slab, p s i, can be evaluated taking into 
account that the thermal pressure in the slab is equal to the magnetic pressure at the 
magneto spheric radius as 

Psi = * T P 6 ■ (16) 
271 Kb To 

Here m p and &b are the proton mass and Boltzmann constant, and Tq is the gas tempera- 
ture at the inner radius of the magnetic slab. 



EVOLUTION OF THE NEUTRON STAR SPIN PERIOD 

The accretion picture in GX 301-2 differs from that expected within the quasi- spherical 
accretion scenario by exceptionally high rate of angular momentum dissipation in the 
accretion flow. We came to this conclusion taking into account that spin-down of the 
neutron star undergoing quasi-spherical accretion is possible only provided the an- 
gular velocity of the accreting material at the magneto spheric boundary, CO em (r m ) = 
££2 rb C^GAm) > is less than angular velocity of the star itself, co s (see Il26ll ). Here 
^orb = 2^/Porb is the average value of angular orbital velocity. Solving inequality 
fflfe m (r m ) < (0 S for we find 

« < 0.03 -- 12/7 ^37 4/7 ^ 0/7 (^) x (17) 

x ( fi VV Vrei Y( B * V /7 

X V685s y \400kms-V \B CRSF J ' 



This value of £ is almost an order of magnitude smaller than the average value of 
this parameter evaluated in numerical calculations of wind-fed accretion made under 
assumption /3 ^> 1 (see [6] and references therein). 

Significant dissipation of the angular momentum in the accretion flow proves to be 
possible if the neutron star is accreting material from the turbulent quasi- statical enve- 
lope rtZZLQj]] or if accreting material possesses strong enough magnetic field [28]. The 



spin-down torque applied to the neutron star from the accretion flow can be estimated as 

= ktMcG^r^, which can be written by substituting M c = /i 2 / (2GM ns r 7 n ) 1/2 in the 
form 

^ (t ) = k t ^i 2 co s = hn 2 
Sd r^ /2 (2GM ns )l/2 (r m r cor ) 3 / 2 ' 



As seen from this expression, the spin-down torque applied to the neutron star from 
the accretion flow essentially depends on its magnetospheric radius and is the greater 
the less the value of magnetospheric radius is. If a neutron star accretes material from 
a quasi-spherical flow or a hot turbulent envelope, its magnetospheric radius is r m ~ r a 



11121. 1271 11 111. However, the spin-down rate of the neutron star in this case turns out to be 
insufficient to explain the observed breaking of GX 301-2 (see expression |5]). 

The magnetospheric radius of the neutron star accreting material from the magnetic 
slab depends on the mechanism of plasma penetration into the stellar magnetic field. If it 
occurs due to interchange instabilities of the magnetospheric boundary (Rayleigh-Tailor 
and Kelvin-Helmholtz instabilities), then, similar to the previous case, the magneto- 
spheric radius will be close to its canonical value, r a . If the interchange instabilities of 
the magnetospheric boundary are suppressed, the plasma penetration into the magnetic 
field occurs due to diffusion process governed by magnetic reconnection at the rate [|29|] : 

M in (r m ) = 4^r m 5 m p Vff (r m ) = 47tr^ 7 l D 1 ^ 'PqV^ 2 \r m ) . (19) 

1 /2 

Here 8 m = (£> e ff Zd) is the depth of the diffusion layer on the magnetospheric bound- 
ary (magnetopause), D e ff is the effective diffusion coefficient and po is the plasma den- 
sity at the boundary. The characteristic time of plasma penetration into the stellar field, 
Zd, is determined by the time on which the plasma having diffused into the field leaves 
the magnetopause moving along the magnetic field lines due to gravitational attraction 
by the neutron star. Under the conditions of interest this time corresponds to the dynam- 
ical (free-fall) time at the magnetospheric boundary, i.e. ~ %(r m ) B29J]. 

It is necessary to note that the situation in which the interchange instabilities of the 
magnetospheric boundary are suppressed due to magnetic shear generation in the mag- 



netopause (see ll30L l31l.l32ll) is not exceptional. The high efficiency of this stabilization 
is well-known, in particular, from laboratory experiments with TOKAMAKs [33]. Be- 
sides, studies of the Earth magnetosphere have shown that the rate at which the solar 
wind plasma penetrates into its magnetic field can be well approximated by the expres- 
sion <fT9l) , in which Z) e ff is the Bohm diffusion coefficient, 

Db = «b^|^, (20) 
2eB(r m ) 

and the efficiency parameter «b ranges between 0.1 —0.25 Q34{] . Here e is the electron 
charge, Tq is the plasma temperature, and B(r m ) is the magnetic field strength at the 
magnetospheric boundary. This provides us with additional grounds to assume that 
plasma penetration into the magnetosphere of the neutron star can occur due to Bohm 
diffusion. 

The magnetospheric radius of the neutron star in this case can be evaluated by taking 
into account the stationary character of the accretion process, that is adopting the rate of 
plasma penetration into the stellar magnetic field equal to the mass capture rate by the 
neutron star from the wind of its companion as well as to the rate of plasma accretion 
onto the neutron star surface , i.e. M- m (r mca ) = LxR n s/GM ns . Solving this equation, we 
find 

r mca * 8 x 10 7 cm x a 2 / 13 ^ R 6 ^\ (21) 



where OTo.i = OCb/0.1 and = To/10 6 K is the plasma temperature at the inner radius 
of the magnetic slab normalized according to the observational results presented in [|35ll . 
The spin-down torque applied to the neutron star from the magnetic slab turns out to be 

(si) _ hn 2 

y mca ' cor ) 

This result makes it possible to estimate an expected rate of the neutron star spin- 
down within the magnetic accretion scenario using the expression V^T = K$ /2%I. 
Substituting the parameters of the pulsar GX 301-2, we get: 

v s ( d mca) ~ TxlO-^Hzs-^ta-f^m- 14 / 13 /^ 1 ^ 13 ^ 3 (23) 



D 57/13 

x R, 



V685s y V 



17/13 



6 V685 S ; V^crsf; 

Thus, the observed spin-down rate of the neutron star in the X-ray pulsar GX 301-2 can 
be explained in terms of the magnetic accretion scenario provided k t > 0.14. 



CONCLUSIONS 

The magnetic accretion scenario has been up to now applied exclusively to the case of 



applied 

accretion on to the black holes (see, e.g. [|13L I17L I18L Il9l \2QX\). It is difficult to perform 



observational checks of main predictions of these studies. Our paper presents one of 
the first attempts to apply the magnetic accretion scenario to reconstruction of the mass 
exchange picture between the components of the X-ray binary system as well as for 
interpretation of observational appearance of these objects. The most important result 
of our study is the conclusion that the spin-down rate of the neutron star undergoing 
magnetic accretion proves to be higher than in case of stars accreting non-magnetized 
material. Observational verification of this result does not encounter any difficulties 
nowadays. This is illustrated by presented above explanation of exceptionally high spin- 
down rate of the neutron star in GX 301-2. This result, however, points out a necessity 
to analyze in more detail the mechanisms used in the magnetic accretion scenario 
which remained poorly understood so far. It is worthwhile to note that the influence 
of the magnetic field of the accretion flow on both its structure and parameters of the 
neutron star magnetosphere can lead to increase of the rate of matter outflow from the 
magneto spheric boundary considered in the paper by Lovelace et al. [360 .The neutron 
star spin-down rate in this case can be higher. Besides, investigating the parameters 
of outflowing matter by means of X-ray spectroscopy opens additional opportunity to 
identify the magnetic accretion scenario. 
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